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Topological gap states of semiconducting armchair graphene ribbons 
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In semiconducting armchair graphene ribbons a tensile strain can induce pairs of topological gap 
states with opposite energies. Near the critical value of the deformation potential these kink and 
antikink states become almost degenerate with zero energy and have a fractional charge of one-half. 
Such a semiconducting armchair ribbon represents a one-dimensional topological insulator with 
nearly zero energy end states. Using data collapse of numerical results we find that the shape of 
the kink displays an anomalous power-law dependence on the width of the local lattice deformation. 
We suggest that these gap states may be probed in optical measurements. However, ’’metallic” 
armchair graphene ribbons with a gap induced by many-electron interactions have no gap states 
and are not topological insulators. 
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FIG. 1: (a) Armchair ribbon with length L y and transverse 
width L x . A and B denote carbon atoms. In the stripe region 
tensile strain that has reflective symmetry about the x axis 
is applied along the ribbon direction, (b) In the stripe region 
the pseudospin density E(y) of a gap state rotates about the 
y axis with only x and 2 components. It represents a kink. 
Outside the stripe region the length of the pseudospin goes 
slowly to zero. 


Introduction - Recently a rapid progress has been made 
in the fabrication of armchair graphene ribbons (TUI] . 
They have a great potential for spintronic applications, 
where electronic nrany-body interactions may play a sig¬ 
nificant role@, 0 ]- Armchair ribbons may also provide 
numerous interesting issues in fundamental physics, such 
as topological objects. By perturbing the graphene lat¬ 
tice and the tunneling coupling constants one can give 
rise to effective gauge and scalar fields p| that generate 
solitonic objects [9jil3J. Generation mechanism of gap 
and vortices is at the center of the study of graphene 
material9, ITol. H34l6| . 

According to tight-binding calculations an armchair 
ribbon has a gap and is semiconducting when the trans¬ 
verse width is L x = (3 M + l)ao or 3Mao, where ao 
is the unit cell length and M is an integer. But when 
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L x = (3 M + 2)ao the armchair ribbon has no gap and 
is metallic. However, the inclusion of many-electron in¬ 
teractions produces a gap@. It is unclear whether these 
armchair ribbons are topological insulators. 

In polyacetylene a fermion mass potential produces 
a gap, and a twist in this scalar potential produces a 
zero energy soliton with mixed chirality and fermion 
fractionalization[9, [HI [l2], [l^]. Similarly, a chiral gauge 
field can simulate the effects of armchair edges and can 
open a gap. However, the Berry phase is absent in arm¬ 
chair ribbons (this can be shown by integrating the tight- 
binding wavefunction over the Brillouin zone). It is thus 
unclear whether zero energy solitons can exist in such 
a case, i.e., whether the armchair ribbon represents a 
one-dimensional topological insulator [l^- The purpose 
of the present paper is to find a distortion of the chiral 
gauge field that will produce a kink, antikink, soliton, 
and charge fractionalization, as the twist of a scalar po¬ 
tential does in polyacetylene. 

Topological properties can be displayed in the variation 
of the pseudospin. Graphene armchair ribbons consist of 
A-carbon and B-carbon atoms, see FigJT| The sublattice 
spin is defined so that when an A-carbon atom is occu¬ 
pied its value is ’up’ and when a B-carbon atom is occu¬ 
pied it is ’down’ (massless fermions can either have their 
sublattice spin pointing along their direction of wave vec¬ 
tor or opposed to it). The valley spin can be defined 
similarly. In an armchair ribbon the sublattice and val¬ 
ley degree of freedom mix, and their spins can be put 
together using 4x4 pseudospin matrices fl8| 

E = /®a, (1) 

where a = ( a x , a y ,a z ) are Pauli spin matrices. We define 
the pseudospin density £;(?/) = 4'*(r)E i 4'(r), where ’l'(r) 
is an eigenstate of the Hamiltonian. A quantum state 
with finite occupation probabilities only on A-carbon 
(’up spin’) atoms to the left of the stripe and only on 
B-carbon (’down spin’) atoms to the right would consti¬ 
tute a topological object called a kink. 

We find that the tensile strain with reflective symme¬ 
try about the x axis can generate a pair of gap states 
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FIG. 2: (a) Eigenenergies for stripe width w = 5A, the num¬ 
ber of fc-space grid points IV = 41, and the strength of de¬ 
formation potential Vf = 2eV (in this paper L x = 25oo and 
L y = 1504A unless stated otherwise. Here ao = 2.46 A is 
the length of the unit cell). Label 1(—1) in the gap denotes 
kink (antikink) with lower (higher) energy, (b) Kink and an¬ 
tikink energies as a function of deformation potential Vf for 
L x = 25ao. Near the critical value Vf they become almost 
degenerate. For larger (smaller) values of L x the energy gap 
is smaller (larger). 


with opposite energies when it is applied along the lon¬ 
gitudinal direction of the ribbon in the stripe region 
—w/2 < y < w/ 2, see Figjl] We find that the properties 
of these states cannot be described entirely by any one 
soliton model (there are several models, namely, the well- 
known free Klein-Gordon equation with a potential with 
a degenerate vacuum[f|, the polyacetylene model[K|, and 
Jackiw-Rebbi solutions at the interface between regions 
of positive and negative masses Q). 

Deformation induced chiral vector potential- In the 
presence of a deformation induced chiral gauge vector 
potential Af the effective mass Hamiltonian of two- 
dimensional Dirac electrons takes the form 

Hyl = v F a ■ (p- -Af(f )) (2) 

for the K valley, and 

H k > = v F cr' ■ (p + -A f (r)) (3) 

for the K' valley (e > 0). Here we have chosen the val¬ 
ley vectors as K = (^,0) and K' = (— ^,0). The 
x and y components of a' are — cr x and a y . Note that, 
unlike for the real vector potential, the chiral vector po¬ 
tential appears with opposite signs [19] for K and K' val¬ 
leys (it is chiral in the sense that it distinguishes val¬ 
leys). Since K and K' valleys are coupled by the arm¬ 
chair edges these two 2x2 Hamiltonians are put to¬ 
gether as a 4 x 4 Hamiltonian, and the resulting Hamil¬ 
tonian is denoted as Hq @,0. For a semiconducting rib¬ 
bon the lowest energy conduction subband has the quan¬ 
tized value k x = -rff— [7] (the metallic ribbon has k x = 0 
and will be treated separately later) The gap 2A sepa¬ 
rating conduction and valence bands can be computed 
from A = ^taok x = 0.098eV (the hopping parameter is 
t = 2.7eV). 


The tensile strain will induce changes of C-C bonds 
along the ribbon direction, in addition to changes of 
the C-C bonds along other directions. In the effective 
mass approximation (ill l 2 h | it produces the chiral vec¬ 
tor potential Af{r) = (AfQ(y), 0,0), where 0(y) = 1 
in the stripe region and 0 outside it, and Af tX = 

£71 — 5(^72 +<^ 73 ) (£71 represents the change of the hop¬ 
ing integral of the bonds along the ribbon direction while 
S 72 and (>73 represent those of side bonds). There is no 
z-component of Af(r), and the y-component is also zero: 
^f-Af^y = f(Ay 2 — < 173 ) = 0 since £72 = £73 for our 
tensile strain. The z-component of the induced magnetic 
field is non-zero only at the edges of the stripe y = ztw/2: 
Bf = Af(S(y + w/2 ) — 8(y — w/2)). This chiral vector 
potential leads to the correction| 22 ] 

V = -V f e(y)Z x ( 4 ) 

in the Hamiltonian. The strength of the deformation is 
Vf = ' A *- F > 0 and has units of energy. Increasing Vf 
corresponds to weakening of the C-C bonds. 

An eigenstate of the total Hamiltonian Hq + V can 
be solved as a linear combination of the basis states 
ip s ,k y if) that are solutions of the Dirac Hamiltonian Hq 
of the armchair ribbon: 'F(r) = J2 s .k y Cs,k y ^ s ,k v (r), 
where C sFy are the expansion coefficients of the lowest 
energy conduction subband (s = 1 ) and highest energy 
valence subband (s = —1). The momentum space inter¬ 
val — k c < k y < k c is divided into N grid points (periodic 
boundary condition is used with k c = -f-N, where N is 

±Jy 

an odd number and L y is the length of the ribbon). The 
results presented below are obtained by this method, or 
by a tight-binding calculation including many-electron 
interactions. 

Kink, antikink, and soliton of semiconducting arm¬ 
chair ribbons.- We show that kink and antikink states 
can exist in the gap. We solve the effective mass Hamil¬ 
tonian Hq in the presence of the deformation potential 
and find that the gap states 1 and \&_ 1 are present 
with opposite energies, see in Figj2] The wave-function 
components of iHi and 4/_i are related through chi¬ 
ral operation: = (T^, fUs, ^' A , 4" S ) T , and T-i = 

(Ta, — 4/s, 4 /( 4 , — ^'b) T ^ where the A- and R-components 
satisfy ^A(y) = —4 , s(— y) and the overall phase factors 
are removed. The probability density of an antikink state 
\F_i is shown in Fig[3Ka). We have computed the angle 6 
between the x axis and S(y), see Fig[3Kb) (The y compo¬ 
nent of the pseudospin density is Ti v (y) = 0). The mag¬ 
nitude of the pseudospin density |E(y)| decays exponen¬ 
tially, see Fig J3](b). In this exponentially decaying region 
the angle is well approximated by 0{y) = ±utanh(ay). 
Here the minus sign corresponds to the antikink state 
while the plus sign corresponds to the kink state. In 
the strong coupling region Vf > 2 eV the amplitude is 
v « 7t/2. 

The energy of a kink is determined by the values of 
v and a, which are related to the physical properties of 
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FIG. 3: (a) Probability density of an antikink state with 
energy E = 0.020eV. The A component (solid line) and B- 
component (dashed line) are shown for N = 81 and Vf = 2eV. 
(b) Numerical result of 9(y) (solid line) for the antikink state 
of (a). Red dashed line is the fitted result using the ampli¬ 
tude v = —1.35 and inverse length a = 0.143A -1 . Also dis¬ 
played is the magnitude of |E(j/)| (dotted line), (c) Compar¬ 
ison between numerical and analytical energies is displayed, 
(d) Data collapse of numerical data. Results are for the stripe 
widths w = 5A (open symbols) and w = 7 A (filled symbols). 
Strengths of deformation potential are Vf = 0.87eV (smaller 
slope) and 2eV (steeper slope). Numbers in the box indicate 
numbers of fc-space grid points N ( k c is proportional to N). 


the kink, namely its amplitude and the inverse size. For 
Vf /Vf < 2 our numerical energies agree approximately 
with the analytical result 

A 8 A 3 

E-i = — - ~(va)£v andFd = + ~(va)fv. (5) 

In Figl3] (c) the energy correction |(ua)^u[§] is shown 
and compared with the numerical result. The parameter 
f = 0.150eVA can be determined by fitting the analyt¬ 
ical energy E to the numerical energy. The deviation 
between the numerical and the analytical energies is the 
largest near the critical value Vf = 2.51eV, which is the 
maximum value of strength of the deformation potential 
Vf « t. Using a tight-binding approach we include the 
effects not included in our model, and find more than two 
gap states. 

Near the maximum value of strength of the deforma¬ 
tion potential Vf two nearly zero energy kink and an¬ 
tikink gap states fki and T_i exist, see Fig©b). We 
will call these kink and antikink solitons. Their proba¬ 
bility densities are similar to the one shown in Fig|3](a). 
These modes are nonchiral. The usual charge counting 
argument gives that the sum of the charges of the two 
gap states is one, implying that each has a fractional 
charge of one halfjUJ. In fact, they are similar (24| to the 
nearly degenerate end states of a one-dimensional topo¬ 
logical insulator with a finite length[25|, despite that the 


Berry phase is zero in the unperturbed semiconducting 
armchair ribbon. 

Anomalous scaling - We find that the field 9(y) satisfies 
near y = 0 the following finite-size scaling j2(| relation as 
a function of the width w and the momentum cutoff k c 


o(y) 

X 0 


= -Ay, 


( 6 ) 


where y = is the dimensionless coordinate and \ = 
is the dimensionless width. From our data collapse, 
shown in Fig©d), and the scaling relation Eq.(© we find 
that the slope of 6(y) displays an anomalous dependence 
on the width w : 


va 


K 2tt J 


(7) 


with the numerical value of (3 = 0.203. The anoma¬ 
lous exponent [27j /? may also be measured experimen¬ 
tally near Vf = 2eV, where the amplitude is v ~ n/2. 
The value of the slope va can be determined from the 
measured values of the kink and antikink energies, given 
by Eq. ©• Then a plot of the slope va as a function of 
the width w would determine the exponent. Also we find 
that the constant A increases as the strength of the de¬ 
formation potential increases. However, as the transverse 
length L x increases, A is approximately unchanged. 

Metallic armchair ribbon.- The inclusion of many- 
electron interactions produces a gap in a metallic arm¬ 
chair ribbon To assess qualitatively the interplay be¬ 
tween the chiral vector potential and many-electron in¬ 
teractions we model it in a half-filled Hubbard model 
of a honeycomb lattice [28|. Applying the unrestricted 
Hartee-Fock method near Vf = Vf, where the C-C bonds 
are cut, we find a significant gap when the on-site inter¬ 
action U is sufficiently large but no gap states (local¬ 
ized states). This is because the on-site interaction U 
shifts zero-energy states into the energy continuum where 
wavefunctions become delocalized. The Kekule-like dis¬ 
tortion, described by 7 °a[I3i, also produces a gap in the 
metallic armchair ribbon. However, adding a chiral vec¬ 
tor potential to it does not generate gap states. The 
metallic armchair ribbon is thus not a topological insula¬ 
tor. In contrast, gap states exist for the semiconducting 
armchair ribbon when U is not too large. 

Discussions and summary - The optical transition from 
a kink H/i to antikink has a substantial opti¬ 

cal strength |(\E , i|j y |'F_i}| 2 when photons are polarized 
along the armchair ribbon ( j y is the current operator). 
It decreases as Vf increases. However, when photons are 
polarized along the transverse direction of the ribbon it is 
zero since |(^i liJlk-i )l 2 = 0. The Coulomb self-energy 
and vertex corrections [29J renormalize the bare transition 
energy upwards about 15%. Optical resonant oscillations 
between H/i and would be also interesting to observe 
near Vf. 

We have investigated gap states of semiconducting 
armchair ribbons within an effective mass Hamiltonian 
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model, supplementing it with a tight-binding calculation 
including many-electron interactions. A chiral deforma¬ 
tion can form kink and antikink pairs with opposite en¬ 
ergies. Although the Dirac singularity at zero energy 
is destroyed by the gap the zero energy solitons, with 
a fractional charge of one half, can exist at the criti¬ 
cal value of the deformation potential. Semiconducting 
armchair ribbons can thus represent a one-dimensional 
topological insulator. Using data collapse of numerical 
results we find that the shape of the kink displays an 
anomalous power-law dependence on the width of the lo¬ 
cal lattice deformation. In contrast ’’metallic” armchair 
ribbons with a many-body induced gap do not have gap 
states and are not topological insulators. We propose 
optical absorption measurements to explore the proper¬ 
ties of a kink, antikink, and solition of semiconducting 


armchair ribbons. Scanning tunneling microscopy might 
also be used to identify the local density of states of the 
topological gap states. Fractionally charged fermions can 
also exhibit non-Abelian braiding statistics [30]. 
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